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ABSTRACT 


DNA barcoding of Bryozoa or “moss animals" has hardly advanced and lacks reference sequences for correct 
species identification. To date only a small number of cytochrome c oxidase subunit I (COI) sequences from 
82 bryozoan species have been deposited in the National Center for Biotechnology Information (NCBI) 
GenBank and Barcode of Life Data Systems (BOLD). We here report COI data from 53 individual samples of 
29 bryozoan species collected from Korean seawater. To our knowledge this is the single largest gathering of 
COI barcode data of bryozoans to date. The average genetic divergence was estimated as 23.3% among 
species of the same genus, 25% among genera of the same family, and 1.7% at intraspecific level with a few 
rare exceptions having a large difference, indicating a possibility of presence of cryptic species. Our data 
show that COI is a very appropriate marker for species identification of bryozoans, but does not provide 
enough phylogenetic information at higher taxonomic ranks. Greater effort involving larger taxon sampling 
for the barcode analyses is needed for bryozoan taxonomy. 
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INTRODUCTION 


Bryozoans or “moss animals" are distributed around the world 
and presently more than 5,000 species are known to exist. 
They are fouling animals, which are easily introduced into 
alien regions via ships and/or ballast waters (Seo, 2005). The 
attachment of these animals to either naturally or artificially 
submerged surfaces, such as ship hulls, enables their frequent 
relocation. Although bryozoans are taxonomically and eco- 
logically important marine invertebrates, DNA barcoding of 
them are not yet progressed. Studies based on the standard 
maker of DNA barcoding of most animal species, mitochon- 
drial cytochrome c oxidase subunit I (COI) sequences, are 
very rare in bryozoans and confined to a few species with, 
for instance, concerns of phylogeography (Gómez et al., 2007) 
and invasion (Mackie et al., 2006). 

Approximately, 130 species are known to be distributed 
in Korea; however, morphological identification of these taxa 
is difficult due to the lack of distinguishing characteristics, 
particularly among closely related species. For this reason, 


we expect that there should be a number of specific synony- 
my. DNA barcoding, therefore, should contribute significant- 
ly to the proper identification of bryozoan species. To date 
DNA barcodes of only 82 bryozoan species have been depo- 
sited in the National Center for Biotechnology Information 
(NCBI) and the database of Barcode of Life Data Systems 
(BOLD). Thus, there is an urgent need for augmenting the 
worldwide database of COI barcodes. 

We attempted to construct COI barcodes based upon 53 
individual samples from 29 species of Bryozoa collected 
from Korean seawater (Table 1). Genomic DNA was isolated 
from a small piece of each species using DNeasy Blood and 
Tissue Kits (QIAGEN Inc., USA) following the manufac- 
turer's protocol. A partial DNA fragment of mitochondrial 
COI was first amplified by PCR using the primers LCO1490 
and HCO2198 (Folmer et al., 1994); however, the PCR pro- 
ducts were not distinguishable on agarose gel. Thus, we de- 
signed a nested second PCR reaction using a degenerative 
BR bugul COI F (5’-CTG GDA TAR TTG GWA GAG 
G-3’) and BR. bugul COI R (5’-GTR TTW AAA TTW 
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Table 1. Classification of Bryozoa (class Gymnolaemata) and GenBank accession numbers 


. Species No. of : 
Order Family genus (GenBank accession no.) individuals Locality 
Cheilostomatida Beaniidae Beania Beania hexaceras 2 Isl. Dokdo 
(HQ896142-3) 
Beania mirabilis 1 Jeju Hwasun Port 
(HQ896144) 
Bugulidae Bugula Bugula neritina 2 Sangju 
(HQ896145-6) 
Bugula subglobosa 1 Isl. Jeju Moon 
(HQ896147) 
Bugula umbelliformis 2 Isl. Goheung Gye 
(HQ896148-9) 
Bicellariella Bicellariella fragilis 2 Isl. Jeju Soop 
(HQ896150-1) 
Candidae Tricellaria Tricellaria occidentalis H Namhae Sangju 
(HQ896152) 
Caberea Caberea lata 1 Tongyeong 
(HQ896153) 
Amastigia Amastigia rudis 2 Isl. Jeju Moon 
(HQ896154-5) 
Amastigia xishaensis 2 Isl. Jeju Soop 
(HQ896156-7) 
Scrupocellaria Scrupocellaria diadema 4 Isl. Jeju Soop 
(HQ896158-61) 
Cellariidae Cellaria Cellaria punctata 2 Isl. Manjae 
(HQ896163-4) 
Hippopodinidae Metroperiella Metroperiella montferrandii 2 Isl. Goheung Doorok 
(HQ896165) 
Metroperiella parviavicularia 1 Isl. Goheung Moohak 
(HQ896166) 
Lanceoporidae Emballotheca Emballotheca pacifica 2 Isl. Jeju Moon 
(HQ896168-9) 
Calyptotheca Calyptotheca symmetrica 2 Pusan Daebyeon Port 
(HQ896170-1) 
Calyptotheca wasinensis 3 Isl. Jeju Moon 
(HQ896172-4) 
Celleporariidae Celleporaria Celleporaria brunnea 3 Pusan Daebyeon Port 
(HQ896175-7) 
Membraniporidae Membranipora Membranipora serrilamella 1 Tongyeong 
(HQ896178) 
Biflustra Biflustra perfragilis 2 Isl. Geoje We 
(HQ896179-80) 
Biflustra crenulata 2 Isl. Geoje Angyeong 
(HQ896181-2) 
Microporellidae Fenestrulina Fenestrulina malusii 2 Pusan Daebyeon Port 
(HQ896183-4) 
Petraliellidae Hippopetraliella Hippopetraliella magna 2 Isl. Somaemul 
(HQ896185-6) 
Exochellidae Escharoides Escharoides excavata 2 Isl. Geoje Jisim 
(HQ896187-8) 
Smittinidae Parasmittina Parasmittina contraria 2 Isl. Goheung Gamnangyeo 
(HQ896189-90) 
Watersiporidae Watersipora Watersipora subtorquata 2 Namhae Sangju 
(HQ896194-5) 
Ctenostomata Flustrellidridae Flustrellidra Flustrellidra armata 1 Isl. Geoje We 
(HQ896196) 
Mimosellidae Bantariella Bantariella bocki 2 Isl. Jeju Moon 
(HQ896197-8) 
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CGR TCK GTT A-3’) primer set specifically designed for 
the order Cheilostomatida, of which COI sequences were 
relatively abundant in GenBank. Fortunately, this new primer 
set was successful for cross-species amplification of COI from 
two species of the order Ctenostomata. The PCR products 
were purified with a LaboPass™ PCR purification kit (Cosmo 
Genetech Inc., Korea) and sequenced using an ABI 3730xl 
(Applied Biosystems Inc., USA). When we found mixed COI 
haplotypes due to a mixture of different individuals of cons- 
pecies, a typical problem of colony forming tissues in some 
specimens, we performed clone sequencing to phase out the 
ambiguous DNA sequence sites. Finally, pairwise sequence 
divergences and a neighbor-joining (NJ) tree were investi- 
gated with a Kimura-2-parameter (K2P) model using the com- 
puter software MEGA 4.0 (Kumar et al., 2004). GenBank 
accession numbers of the sequences that we deposited from 
this study and their taxonomical information are provided in 
Table 1. 


RESULTS 


Our COI data was built from the largest taxon sampling to 
date, representing 29 species, 18 genera, 15 families, and 2 
orders of Bryozoa. Assuming that the present COI data were 
obtained from the recognized species without error, the esti- 
mates of pairwise sequence divergence between species are 
as follows: 25% at a confamilial level, 23.3% at a congeneric 
taxonomic level, and 1.7% at a conspecific level (Fig. 1). 
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Fig. 1. Distribution of pairwise sequence divergence based on 
the Kimura-2-parameter for cytochrome c oxidase subunit I 
sequences according to taxonomic levels. The horizontal axis 
represents intervals of genetic distance in percentage and the 
vertical axis is the number of individuals associated with each 
distance interval. 
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Bryozoan DNA Barcode 


Although Fig. 1 highlights increasing genetic divergence as 
taxonomic rank increases, the difference between generic 
and familial levels is still very small. Therefore, the main 
purpose of DNA barcoding for species identification has been 
fulfilled given the current taxon sampling with some excep- 
tionally high divergent conspecific cases which strongly in- 
dicates the possibility of the presence of cryptic species in 
Scrupocellaria diadema, Watersipora subtorquata and Para- 
smittina contraria. Further examinations are needed for these 
cases. The NJ tree clearly shows high bootstrapping supports 
in most terminal nodes, which indicates that COI sequences 
have sufficient information for the identification of species 
given the present taxon sampling. Among all 29 species, only 
two, Biflustra perfragilis and Biflustra crenulata, could not 
be distinguished due to identical sequences. In contrast, we 
could not find such high bootstrap supports in most basal 
nodes of the tree, suggesting a lack of phylogenetic signal at 
higher taxonomic levels (Fig. 2). Furthermore, most families 
with multiple genera were not grouped as monophyletic 
clades, and many basal branches show polytomies again indi- 
cating a lack of phylogenetic signal. 


DISCUSSION 


Our study of bryozoans demonstrates the classic problems of 
species identification by molecular markers, which suffers 
from insufficient samplings, superficial knowledge of the 
range of molecular divergence, cryptic species, and pheno- 
typic plasticity of recognized species (Vrijenhoek, 2009). In 
some cases, the COI marker did not provide resolution of 
species identification. As an extreme example, four sequences 
of Biflustra perfragilis and Biflustra crenulata turned out to 
be identical, while they are easily classified by distinctive 
morphology. Therefore, it will be interesting to prove whether 
the taxonomic key characteristics for the two species result 
from phenotypic plasticity or the COI marker is just not ap- 
propriate for differentiate such closely related species. On 
the other hand, Scrupocellaria diadema showed a relatively 
large difference among different individuals (n=4) (Fig. 2). 
This species is known to have relatively large number of com- 
plex key characteristics and show a particularly high degree 
of morphological variability (Seo, 2005). Therefore, it is pos- 
sible that this species may be a species complex. Further 
molecular examinations based on multiple individuals and 
localities are clearly warranted for this species. In addition, 
extremely divergent individual sequences were observed 
from Watersipora subtorquata and Parasmittina contraria, 
respectively. Interestingly the genus, Watersipora, includes 
cosmopolitan species and has been the subject of the studies 
of alien species (Seo, 1999). Because this taxon has a rather 
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Fig. 2. Neighbor-joining (NJ) tree of Bryozoa based on cytochrome c oxidase subunit I sequences from 54 bryozoan samples. The 
root was drawn using the midroot method and the asterisks (x) indicate two reference sequences (Securiflustra securifrons, FJ196090; 


and Flustra foliacea, FJ196080) of the order Cheilostomatida 


obtained from GenBank for within group comparison. Numbers on 


branches show bootstrap values (» 5096) in NJ analysis with 500 replicates. 


simple morphology, morphological species identification is 
difficult and, accordingly, multiple cryptic species could be 
involved. In contrast, Parasmittina contraria is known to 
highly variable species, but our COI data strongly suggest a 
possible presence of cryptic species within this taxon. This 
species is currently known to distribute only in the Jejudo 
Islands. The high divergence of Watersipora subtorquata 
and Parasmittina contraria samples corresponds to the right 
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tail of the distribution of pairwise sequence divergence at the 
conspecific level (Fig. 1). Our DNA barcoding clearly sepa- 
rate the morphologically distinctive species including Calyp- 
totheca wasinensis, Emballotheca pacifica, Celleporaria 
brunnea, Bugula umbelliformis, Cellaria punctata, and Ban- 
tariella bocki (Fig. 2). One of the advantages of COI sequ- 
ences and their comparison could be added by the case of 
Fenestrulina malusii. Particularly this species causes a con- 
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siderable difficulty in taxonomic classification due to its am- 
biguous morphological characters. Thus the resultant diver- 
gent COI sequence of it to other recognized bryozoan species 
can be used as one of many references to identify it in future 
(Fig. 2). 

It should be noted that one of the most fundamental pro- 
blems encountered with DNA barcoding of bryozoans was 
the lack of a universal COI primer set and insufficient refer- 
ence sequences. During our study, we experienced frequent 
PCR failures with the universal Folmer primers (Folmer et 
al., 1994) for many bryozoans (not shown in the present re- 
sults). This problem is not surprising if we take into account 
the fact that taxonomically broad bryozoans may have an 
enormous degree of COI divergence. In fact, the range of 
divergence and polymorphism of COI in bryozoans is yet 
unknown. The previously known COI sequences are limited 
to a very small portion of bryozoans, which actually impedes 
the design of universal primers. Hopefully, as COI data of 
bryozoans grow, development of universal primers should 
be possible. Such group specific oligonucleotide sequences 
might be desirable to minimize contamination due to non- 
bryozoan amplification, known as "the peril of universal 
primers" (Vrijenhoek, 2009). Another experimental difficulty 
was the frequent observation of mixed sequence reads from 
well-amplified PCR products, which indicates two possibil- 
ities: that is, either a mixture of variable haplotypes originat- 
ing from a tissue specimen consisted of numerous tiny colo- 
nial individuals or a contamination of quite different species. 
Although the characteristic colonial structure of bryozoans 
imposes technical difficulties for extracting clean COI sequ- 
ences from raw specimens. If we have suitable COI data from 
taxa across broad bryozoans taxa, the problem of mixed sequ- 
ence reads will be more easily solved by comparing reference 
sequences well annotated with taxonomic information. In 
this regard, our current data, except for some suspicious cases, 
will add better discriminating power to the COI database in 
GenBank and BOLD for comparative analyses. 

In summary, our DNA barcode study involving 29 bryo- 
zoan species is the first attempt of its kind in Korea. We found 
a high degree of COI sequence divergence among most spec- 
ies. The amount of differences was sufficient for species 
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identification in most cases except for two species. We there- 
fore believe that COI can be served as a standard molecular 
maker for DNA barcoding of bryozoans. 
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